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ABSTRACT 

Recent high-resolution calculations of neutron star coalescences are used to investi- 
gate whether vV annihilation can provide sufficient energy to power gamma-ray bursts, 
especially those belonging to the short duration category. Late time slices of the sim- 
ulations, where the neutrino emission has reached a maximum, stationary level are 
analyzed to address this question. We find that vV annihilation can provide the neces- 
sary driving stresses that lead to relativistic jet expansion. Maximum Lorentz factors 
of the order of 15 are found along the binary rotation axis, but larger values are 
expected to arise from higher numerical resolution. Yet the accompanying neutrino- 
driven wind must be absent from the axis when the burst occurs or prohibitive baryon 
loading may occur. We argue that even under the most favorable conditions, vV an- 
nihilation is unlikely to power a burst in excess of ^ 10*^ erg. Unless the emission is 
beamed into less than one percent of the solid angle, which we argue is improbable if 
it is coUimated by gas pressure, this may fail to satisfy the apparent isotropic energies 
inferred at cosmological distances. This mechanism may nonetheless drive a precursor 
fireball, thereby evacuating a cavity into which a later magnetic driven jet could ex- 
pand. A large range of time delays between the merger and the black hole formation 
are to be expected. If the magnetic driven jet occurs after the black hole has formed, 
a time span as long as weeks could pass between the neutrino powered precursor and 
the magnetic driven GRB. 

Key words: dense matter; hydrodynamics; neutrinos; gamma rays: bursts; stars: 
neutron; methods: numerical 



1 INTRODUCTION 

Double neutron star (NS) binaries, such as the famous 
PSR1913-I-16, will eventually coalesce due to angular mo- 
mentum and energy losses to gravitational radiation, pro- 
vided that their orbital separation is small enough. When 
a NS binary coalesces, the rapidly-spinning merged system 
could be too massive to form a single NS; on the other 
hand, the total angular momentum is probably too large 
to be swallowed immediately by a black hole (BH). The ex- 
pected outcome would then be a spinning BH, orbited by a 
torus of NS debris, whose accretion can release sufficient 
gravitational energy, « 10^^ erg, to power a gamma-ray 
burst (GRB; Lattimer & Schramm 1976; Paczyhski 1986, 
1991; Eichler at al. 1989; Narayan, Paczyhski & Piran 1992; 
Mochkovitch et al. 1993; Kluzniak & Lee 1998; Rees 1999; 
Rosswog et al. 1999; RufFert & Janka 1999; Salmonson, Wil- 
son & Mathews 2001; Rosswog & Davies 2002). The energy 
released in the accretion process is expected to be trans- 
formed into a fireball, which may provide the driving stress 
necessary for refativistic expansion (see Piran 1999 and 
Meszaros 2002 for recent reviews) . Possibie forms of this out- 
flow include the following. Relativistic particles generated by 



uV annihilation - the thermal energy released as neutrinos is 
reconverted, via collisions, into pairs and photons. Elec- 
tromagnetic Poynting flux - as in pulsars, strong magnetic 
fields anchored in the dense matter convert the rotational 
energy of the system into a Poynting-dominated outflow. In 
either case, the duration of the burst is determined by the 
viscous timescale of the accreting gas, which is substantially 
longer than the dynamical or orbital timescale, thus provid- 
ing a simple interpretation for the large difference between 
the durations of bursts and their fast variability. 
These binary encounters are thought to transpire outside 
star forming regions (Fryer, Woosley & Hartmann 1999) 
and are not currently thought to be appropriate for long 
GRBs (> 20 s), whose afterglows have been predominantly 
localised within the optical image of the host galaxy (Bloom, 
Kulkarni & Djorgovski 2001). However, these calculations 
are uncertain because they are sensitive to a number of pa- 
rameters that are not well understood (e.g. distribution of 
initial separations; see Belczynski et al. 2002). Such mod- 
els may also be in difficulties producing adequate duration, 
coUimation, and energy (Katz & Canel 1996; RufFert et al. 
1997; Popham, Woosley & Fryer 1999; Ruffert & Janka 1999; 
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Lee & Ramirez-Ruiz 2002), but each of these problems is 
model-dependent and might be cured by different assump- 
tions regarding the poorly known input physics such as disk 
viscosity or the jet production mechanism. 
In this Letter, we study the neutrino emission produced in 
the shocked and dissipatively heated merger remnant of a 
NS binary encounter. Late time slices of three dimensional 
(3D), high resolution calculations of the last stages of the co- 
alescence are used (Rosswog & Davies 2002; hereafter RD) 
to investigate the potential of w annihilation as a viable 
source for GRB production, at least for bursts belonging 
to the short category. Neutrinos could give rise to a rela- 
tivistic pair dominated outflow if they annihilate in a 
region of low enough baryon density. To that effect, high res- 
olution simulations seem to be required, since even a very 
small amount of baryons polluting the outflow could severely 
limit the attainable Lorentz factor. Information regarding 
how much i/V energy and momentum has been injected, its 
confinement and coUimation are at the forefront of our atten- 
tion. We address all three issues here, along with the type 
of predictions that would help us to discriminate between 
this mode of energy extraction and that of strong magnetic 
fields anchored in the dense matter. 



2 NUMERICAL METHODS 

The 3D hydrodynamics simulations of the mergers of bi- 
nary NS were performed by RD using a smoothed particle 
hydrodynamics method with up to ~ 10^ particles. A re- 
alistic equation of state for hot, dense nuclear matter has 
been used, which is based on the tables provided by (Shen 
et al. 1998a, b) and smoothly extended to the low density 
regime with a gas consisting of neutrons, alpha particles, 
photons and pairs (see RD for details). This allows 
one to closely follow the thermodynamic evolution of the 
neutron-rich debris after the merger. Under the conditions 
encountered in the debris - temperatures of several MeV 
and densities around ~ 10^^ g cm~'^ - neutrinos are emit- 
ted copiously. Their effect on the cooling and the changes in 
the composition of the material are taken into account via a 
detailed multi-flavour neutrino treatment which is described 
in detail in Rosswog et al. (2002). The results presented here 
are based on the close analysis of late time segments, typi- 
cally a time tsim ~ 15 ms after the start of the merger sim- 
ulations, where the neutrino luminosities have reached their 
maximum, stationary level (Rosswog et al. 2002). Since we 
are interested in the maximum possible effect, we assume 
the system to maintain this level of emission for as long as 
1 s (see §3). We describe results from three representative 
runs. First, an initially corotating system with twice 1.4 M0 
(run D from RD, referred to in the following as cl.4), which 
yields the lowest temperatures of all runs; it can therefore 
be considered as a lower limit on the neutrino luminosities. 
Second, a system with twice 1.4 M0 and no initial NS spins 
( "irrotational" ; run C from RD, hereafter il.4), which we 
regard as the generic case and finally, an extreme system of 
twice 2.0 M© and no initial spins (run E from RD, hereafter 
i2.0), which we consider to be an upper limit for the neu- 
trino luminosity. 
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Figure 1. Colour-coded is the ratio of energy deposited via uV — > 
e+e~ to rest mass energy, r], which is a measure of the maximum 
attainable Lorentz factor. Shown are the values of log(»7) in the x- 
z-plane above the merged remnant of model il.4 (no initial spins, 
2 X 1.4 M0) at the end of the simulation (t=13.8 ms). Due to 
the symmetry of the merger remnant with respect to the orbital 
plane a similar jet will occur along the negative z-axis. 



3 e± FIREBALLS PRODUCED BY uV 
ANNIHILATION 

The vV e^e~ process (which scales with the square of the 
neutrino luminosity) can tap the thermal energy of the hot 
debris disk. For this mechanism to be efficient, the neutri- 
nos must both escape before being advected and not being 
produced too gradually. Typical estimates for coalescing NS 
suggest a limit of < 10^^ erg for the neutrino energy dumped 
into -pairs (Ruffert et al. 1997; Ruffert & Janka 1999; 
Popham et al. 1999) . If the -dominated plasma were col- 
limated into a solid angle O.^-, then of course the apparent 
"isotropized" energy would be larger by a factor of An/Q,^-. 

The presence of a region of very low density along the 
rotation axis, away from the equatorial plane of the de- 
bris must be present when the burst occurs or otherwise 
prohibitive baryon loading will occur - for instance an en- 
ergy of 10^^ erg deposited into pairs could not accel- 
erate an outflow to F > 100 if it had to drag more than 
E/{Vc^) « 5 • IO-^Mq baryons with it. In Fi gure 1 we plot 
the ratio of energy deposition to baryon rest mass en- 
ergy, r) — (5^-Tinj/(pc^), in the region above the poles of 
the merged remnant, p denotes the matter density and, for 
simplicity, an energy injection time, Tinj = 1 s, well within 
the distribution of the short-duration bursts has been as- 
sumed. T) can be understood as an indicative of the terminal 
bulk Lorentz factor, F. The above t) estimates assume, also 
for simplicity, that the energy deposition rate by vV anni- 
hilation into pairs at tsim is both representative of the 
subsequent phases and steady for at least one second (this is 
clearly an upper limit to the total injected energy over the 
assumed period Tinj). Due to the flnite resolution of the sim- 
ulations the densities along the rotation axis are determined 
by particles located in the inner parts of the disk. We there- 
fore expect to overestimate the densities along the rotation 
axis. Given the above assumptions, higher numerical reso- 
lution would yield higher values for r). This is supported by 
calculating these densities for two equivalent runs of differ- 
ent numerical resolution (run A and B of RD). The densities 
in well-resolved case are lower by approximately one order 
of magnitude. 

The resulting constraints on the injected energy as a 
function of r) are illustrated in Figure 2. It shows the fast 
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Figure 2. Energy deposited via i^F — > e+e~ as a function of rj 
for NS-NS merger models cl.4 (corotation, 2 x 1.4 Mq), il.4 (no 
spins, 2 X 1.4 Mq), i2.0 (no spins, 2 X 2.0 M0). 



T 




2 



Figure 3. Volume of the compact region containing material with 
a given rj for NS-NS merger models cl.4,il.4, 12.0. 

drop of the contained energy as a function of the maximally 
reachable Lorentz factor. The copiously emitted electron 
neutrinos and anti-neutrinos dominate the annihilation pro- 
cess by ~ 95 %, their mean energies (E) arc ^ 10 MeV for 
and ~ 15 MeV 17^. The highest annihilation rates occur 
in the uninteresting (since extremely baryon-rich) hot inner 
disk regions. Even in the most optimistic case (i2.0) only 
~ 10** (47r/f2^77) erg arc found able to blow out material 
near the axis with a Lorentz factor larger than 10 (indeed 
no SPH-particles are enclosed in this region, but wo arc lim- 
ited in our density calculation by the finite resolution). A 
broad spread in the Lorentz factor will therefore be present 
in the outflow close to the rotation axis F may bo high 
(i.e. > 10); at large angles away from the axis an increasing 
degree of cntrainmcnt will correspond to a decrease in F (see 
Fig. 1). Even if the outflow is not narrowly collimatcd, some 
beaming is expected because the disk geometry (see Fig. 15 
RD) channels the energy preferentially along the rotation 
axis. But unless is < 10"^ — 10"'^ this may fail to sat- 
isfy the apparent isotropized energies of ~ 10®^ erg implied 
for "short-hard" bursts at a redshift z = 1 (Panaitescu, Ku- 



mar & Narayan 2001; Lazzati, Ramirez-Ruiz & Ghisellini 
2001). 



3.1 Beaming and Confinement 

When a large amount of energy is suddenly released into a 
compact region, an opaque "fireball" is created due to the 
prolific creation of pairs (Cavallo Sz Rees 1978; Goodman 
1986; Paczyriski 1986; Shomi & Piran 1990). A diagram sim- 
ilar to Figure 2 can be drawn to illustrate the volume that 
contains material above a given r] (see Fig. 3). Consider a 
homogeneous fireball of energy S^u, total mass Mq initially 
confined to a compact region, whose volume is given by Vcav. 
Clearly, since the optical depth r > 1, the initial fireball will 
be an opaque sphere in thermal equilibrium, characterised 
by a single temperature: 

« lOS^/^yJ/t MeV, (1) 

where we adopt the convention Q = 10^ Qa;, using cgs units. 
This radiative sphere expands and cools rapidly until the 
energy of the photons degrades below the production 
threshold. When some amount of baryonic matter is mixed 
with the fireball and the radiation energy dominates the evo- 
lution (i.e. rj ^ 1), the fluid expands under its own pressure 
such that its Lorentz factor grows initially linearly with ra- 
dius, F cx r, and reaches a final value of F ~ r; (Piran 1999). 

If the flow is assumed to be stationary and isentropic, 
then the fluid velocity will increase as the external pres- 
sure of the surrounding medium decreases. When the pres- 
sure has halved the flow becomes transonic and the cross- 
sectional area is minimised. In this way a directed de Laval 
nozzle can be established. Unlike the static situation envis- 
aged by Blandford & Rees (1974), the injected in this 
plasma do already have a relativistic outward motion and 
the flow is then unlikely to be recollimated in a fluid nozzle. 
Although it will become free and supersonic (indeed rela- 
tivistic, if it was so initially) since the pressure drop at fcav is 
very steep. In the simulations, the external pressure is found 
to vary with distance along the rotation axis as poxt oc r"~'^. 
If the pressure drops faster than r^^ onwards, sound waves 
will eventually be unable to cross the jet, which will become 
overpressured with respect to the external medium, and 
thenceforth expand out over an angle F"'^. If the free expan- 
sion starts just outside rcav, where F ~ 2, then it will spread 
over a wide angle (i.e. 6^-^ ~ 30°) and develop into a roughly 
semi-spherical blast wave. The flreball will then accelerate 
until the entire energy is converted into kinetic energy at 
r;r cav li t] > rit, ^ l^OB^^^^gr^^^g , the flreball continues to 
be radiation dominated when it becomes optically thin and 
most of the energy escapes as photons. An observer will de- 
tect them with a characteristic thermal peak frequency of 
Tobs ~ T^T7 (Goodman 1986), which is likely to be outside 
the BATSE [20-600] keV spectral window (see equation 1). 
Alternatively, if 77 < 77b, the flreball becomes matter domi- 
nated, and most of the initial energy is converted into bulk 
kinetic energy with F m r] < rjb. The inertia of the swept-up 
external matter will then decelerate these ejecta significantly 
by the time it reaches rr ~ W^^ E^-'^.„ri7^^^6~-{'^n~^^n cm 

iyu,4S '-L uu.l ism,[) 

(Meszaros & Rees 1997). For an approximately smooth dis- 
tribution of external matter, the bulk Lorentz factor of the 
fireball decreases as an inverse power of the time (oc t~^^^). 
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As a consequence, the spectrum softens in time as the syn- 
chrotron peak corresponding to the minimum Lorentz fac- 
tor and magnetic field decreases, leading to the possibility 
of multi-wavelength afterglow radiation. 

3.2 Accompanying winds 

The neutrinos that are omitted from the inner shock-heated 
regions of the torus will deposit part of their energy in the 
outer parts of the thick disk that has formed around the 

central object (sec Fig. 15 of RD). The possibility of such 
a wind from a NS coalescence has already been realized by 
Ruffert et al. (1997). Using the typical numbers from our 
simulations, we find that the neutrinos will drive a mass 
outflow at a rate (Qian & Woosley 1996): 

*...0-"M,/.4£(3Sv)*(^)'(^)^P' 

El 

where Lr^^ is the luminosity in i^e-neutrinos, e^^ = {-^), 
Rem and Mem arc the radius of and the mass inside the 
neutrino emitting region. This estimate is consistent with 
those of other authors (Duncan, Shapiro & Wassermann 
1986, Woosley 1993). This wind accelerates the blown-off 
material to its asymptotic velocities of ~ 10^ cm . A 
wind of this strength represents a possible danger for the 
emergence of a highly relativistic jet. The jet can only ac- 
celerate to high Lorentz factors if the wind is kept out of the 
funnel region along the rotation axis via centrifugal forces. 
However, from the bulged geometry (see Fig. 15 RD) and 
the assumption that the neutrinos will be emitted preferen- 
tially in the direction h = — Vp/| — Vp|, we conclude that 
most of the wind will be still directed away from the jet. 



4 DISCUSSION 

4.1 uV annihilation as a viable energy source? 

We now want to address the viability of NS-NS binary coa- 
lescences as central engines of GRBs. The main form of en- 
ergy release from the disk we considered is neutrino emission, 
which produces relativistic jets but fails to produce very en- 
ergetic bursts. Due to our post-processing approach the in- 
teraction between the jet and the matter evolution/neutrino 
emission is neglected. Not accounting for the evolution of 
the jet during the energy injection period leads to an over- 
estimation of the energy density in the jet. But even under 
the most favourable conditions (maximum baryon free vol- 
ume allowed by the SPH-particle distribution, neglecting a 
neutrino driven wind that might pollute the funnel region 
above the poles, see §3, and ignoring the jet evolution) we 
find typical energies of only ~ 10*^ erg, where we have as- 
sumed a steady injection lasting for as long as 1 s. Unless 
the plasma was coUimated into a solid angle Q^- < 10^^ 
(which is unlikely if the plasma is solely confined by gas 
pressure; but see Lcvinson & Eichler 2000 for an alternative 
mechanism for liydrodynamic collimation) , vV annihilation 
could only be responsible for a rather weak energy release 
(which is also likely to be very short, almost impulsive; see 
Lee & Ramirez-Ruiz 2002) and may therefore fail to pro- 
duce isotropized energies of ~ 10^^ ergs that are estimated 



for "short-hard" bursts if they are at z ~ 1. Even if this is the 
case, the deposited by vV annihilation may be responsi- 
ble for either precursor emission (thermal and non-thermal) 
and/or for creating a cavity into which a magnetised jet can 
subsequently expand (see below). In this latter case, the af- 
terglow, at legist over a range of directions, would not arise 
until the late magnetised ejecta hits the wall of the cavity. 
One would naturally expect a wide variety of afterglow be- 
haviours arising from the impact of (possibly anisotropic) 
ejecta on an irregularly shaped cavity. 
It is, of course also possible to produce more energetic bursts 
via vV annihilation if one is able to increase the temper- 
atures in the appropriate density regimes (i.e. ^ 10^^ to 
10^'^ gcm~^). One such possibility would be that the "real" 
equation of state is much softer than the one used here, 
leading to higher temperatures and hence more energetic 
neutrino emission. The neutrino luminosities are sensitive 
functions of the temperature of the emitting region, the (3- 
reactions are oc and the plasmon decay is (x T^. Higher 
temperatures would translate into higher neutrino energies 
E^, which in turn render the debris more opaque to the 
neutrinos (the cross-sections scaling oc E^). We still, how- 
ever, expect higher luminosities and hence higher annihila- 
tion energies if the temperature is increased. An alternative 
possibility may arise from general relativistic effects, which 
are not included in the current 3D hydrodynamics simular 
tion. The tendency of general relativistic effects would be to 
compactify the matter distribution, and so higher temper- 
atures would be naturally expected. However, in this case 
the annihilation process would be additionally complicated 
by effects such as bending of neutrino trajectories and red- 
shifts of the neutrino energies. These effects have been found 
to partly cancel out one another so that they do not al- 
ter the results substantially. For a further discussion of this 
point we refer to the literature (Salmonson & Wilson 2001, 
Asano & Fukuyama 2000, Jaroszynski 1993). A third pos- 
sibility could arise from a very high disk viscosity. In our 
current simulation no explicit, physical viscosity has been 
introduced. We measured an equivalent a-viscosity of the 
order of adiak « 4 • 10~^ (lower values are found in the bet- 
ter resolved central object; see RD). The reason for this low 
value is a largely improved artificial viscosity scheme and 
the high resolution. Future investigations including explicit, 
physical viscosity should claxify this point. The raised possi- 
bilities should certainly be kept in mind, but at the present 
stage they are a matter of speculation. 

4.2 Delayed BH formation and MHD jets 

An alternative way to tap the debris energy is via mag- 
netic fields threading the torus. Even before a BH forms, 
the merging system might lead to winding up the fields and 
dissipation before the merger. While the above mechanism 
can tap the rotational energy available in the debris torus, a 
BH formed from the coalescence is guaranteed to be rapidly 
rotating, and being more massive, could contain a larger en- 
ergy reservoir than the debris itself. A magnetic configura- 
tion capable of powering bursts requires fields of a few times 
10^'"" G. A weaker field would extract inadequate power; how- 
ever, it only takes of order of a second for simple winding 
to amplify a 10^^ G field to 10^^ G (Kluzniak & Ruderman 
1998). If magnetic fields of comparable strength than those 
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anchored in the torus thread the BH, its rotational energy 
offers an extra (an possibly even greater) source of energy 
than in principle can be extracted via the B-Z mechanism 
(Blandford & Znajek 1977). For a maximally rotating BH, 
this is 0.29 Mbh(? erg, multiplied by some efficiency factor. 
Even allowing for low total efficiency (~ 20%), a jet powered 
by the B-Z effect would not require any beaming to produce 
the equivalent of an isotropic energy of 10^^'^ erg. The en- 
trained baryonic mass would then only need to be below 
1O~''M0 (which is clearly the case along the rotation axis 
of the merger remnant) to allow high relativistic expansion 
speeds. 

The masses of the central objects of the merger remnant 
are ~ 2.5 and therefore above the value of ~ 2.2 M© 
that is typically found for cold, non-rotating NS in beta- 
equilibrium (see e.g. Akmal et al. 1998). Due to the poorly 
known physics under these extreme conditions (strong field 
gravity, magnetic fields, 'exotic' matter in the NS cores) the 
time scale on which the central object is going to collapse 



is, how(tvpr uncertain For the generic NS merger case^ i1 4 
there are no signs of a collapse - at the end of the simulation, 
the maximum density of the merger remnant has not even 
reached the central density of a cold, non-rotating single 1.4 
Mq NS. This is mainly due to the stabilising effect of differ- 
ential rotation. For a discussion of further stabilising effects 
(some of which are not included in the simulation) we refer 
to RD. The merger remnant would lose angular momentum 
on a magnetic dipole radiation time scale 



10^ s 



Mco 
2.5 Mp 



10^" G 
B 



which is substantially shorter than the viscous timescale, 10^ 
s (Shapiro 2000). Due to the dependence of this timescale 
on the input quantities, we expect a large range of possible 



delays jetween the merger and the subsequent BH forma- 
tion. This will be mainly determined by the details of the 
coalescence, unless exotic matter appears already for very 
low masses, therefore leading to an immediate collapse. Note 
that with reasonable assumptions even a time scale of weeks 
is still possible. More prominent time delays could also be 
expected for mergers of unequal mass NS, or NS with other 
compact companions. Such events might also lead to repeat- 
ing episodes of accretion, or to the eventual explosion of the 
NS which has dropped below the critical mass, all of which 
would provide a longer lasting energy output. A GRB pro- 
duced by MHD coupling is likely to go off inside a cavity 
inflated by the vV annihilation burst (which may not pro- 
duce a standard GRB; see §3). Such cavities can be as large 



as fracliiuus uf a paisyc, giving lisy Lu a maguyLic decelyia- 
tion sh ock montha after the firat vv annihilation burat. 



Much progress has been made in understanding how gamma- 
rays can arise in fireballs produced by brief events deposit- 
ing large amounts of energy in a small region, and in deriv- 
ing the generic properties of the corresponding afterglows. 
Nonetheless, there still remain a number of mysteries, espe- 
cially concerning the progenitors of short-duration bursts. 
As argued above, a wide diversity of behaviours may be the 
rule, rather than the exception. And if the class of models 
that we have advocated here turns out to be irrelevant, the 
trigger of these events has to be even more remarkable and 
exciting. 
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